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Abstract 

We study a generalized two-species model on a ring. The original model [1] 
describes ordinary particles hopping exclusively in one direction in the presence 
of an impurity. The impurity hops with a rate different from that of ordinary 
particles and can be overtaken by them. Here we let the ordinary particles hop 
also backward with rate q. Using the Matrix Product Ansatz (MPA), we obtain 
the relevant quadratic algebra. A finite dimensional representation of this algebra 
enables us to compute the stationary bulk density of the ordinary particles, as well 
as the speed of impurity on a set of special surfaces of the parameter space. We 
will obtain the phase structure of this model in the accessible region and show how 
the phase structure of the original model is modified. In the infinite-volume limit 
this model presents a shock in one of its phases. 
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1 Introduction 



Recently much attention has been focused on one-dimensional reaction-diffusion pro- 
cesses. These models can describe many physical phenomena such as hopping conductiv- 
ity, growth processes and traffic flows [2-4]. They are also of interest from the mathemat- 
ical point of view due to their relation to integrable quantum chain Hamiltonians [5,6]. 
The simplest model of this kind is the Asymmetric Simple Exclusion Process (ASEP) 
with open boundaries [7]. This model comprise particles which jump independently to 
their right with hard core repulsion along a one-dimensional lattice. The open boundary 
conditions mean that particles are injected at one end of the lattice and arc removed at 
the opposite end. This model exhibits a shock structure in the density profile of parti- 
cles. In periodic boundary conditions case the microscopic location of the shock can be 
identified by defining a second class particle (impurity). The ASEP in the presence of an 
impurity on a ring has been studied in the two following cases: 

1) In the first case, the single impurity hops in the opposite direction relative to the 
ordinary particles [8,9]. In this case a first-order phase transition between a low-density 
and a traffic jam phase is observed . 

2) In the second case, the impurity moves in the same direction as the ordinary par- 
ticles [1,17]. The phase diagram of this model consists of six distinctive phases (I- VI) 
in which two of them are symmetric to other phases under a charge conjugation and 
reflection symmetry. The authors have shown that one phase exists in the system in 
which the impurity causes a shock. 

Another example of such driven diffusive systems is the Partially Asymmetric Simple 
Exclusion Process (PASEP). In this model the particles are allowded to jump both to 
their immediate right (with rate p) or left (with rate q) site, if the target site is not 
already occupied. This model has been extensively studied both with open boundaries 
and on a ring [10]. 

In this paper we will study the effects of the presence of a single impurity on the 
PASEP on a ring. Here the ordinary particles can hop to their immediate right (left) 
site, provided that it is empty, with rate 1 (g < 1). The single impurity can only hop 
to its immediate right site (if it is not already occupied) with rate a{< 1) and can 
be exchanged from the left with the ordinary particles with rate /3(< 1). For q = 
this model reduces to the model (2) as discussed above. Using Matrix Product Ansatz 
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(MPA) introduced in [11], we obtain the relevant quadratic algebra which has both finite 
and infinite dimensional representations. For simplicity, we adopt a finite dimensional 
representation of the algebra and carry out all the calculations using a grand canonical 
ensemble in which the population of ordinary particles can fluctuate. By adjusting the 
fugacities of the ordinary particles and of the holes, one can produce some fixed densities 
for them. Although the finite dimensional representation restricts us to the region under 
the surface a -\- 13 -\- q — lin the three-dimensional parameter space, nevertheles we 
shall find the exact phase structure and calculate precisely the density profile of ordinary 
particles and the speed of impurity in this region. With this exact results, we will show 
that three phases exist in this region. In two of them, which are symmetric to each other 
under a charge conjugation and reflection symmetry, the density profile of particles has an 
exponential behavior. We will determine the relevant correlation lengths in these phases 
and the critical values of the rates that characterize the divergence of these correlation 
lengths. In the third phase the density profile of the ordinary particles is linear which is 
a signature of a shock. 

This paper is organized as follows. In section 2 we will describe the model. In section 
3 we will write the weights of the configurations in the stationary state in terms of a trace 
of L + 1 non-commuting matrices and obtain the quadratic algebra. In section 4 we will 
introduce all possible representations of the quadratic algebra and discuss different phases 
in the accessible region of the phase space. In section 5 using one and two-dimensional 
representations of the algebra we will obtain some exact results for an infinite system. 
In the last section we will compare our results with those obtained in [1] for g = 0. 

2 The Model 

The model proposed here contains two species of particles on a closed ring of L + 1 sites 
which are labled from 1 to L + 1. We specify each configuration of the system by an 
L + 1-tuple (ri,r2, t/^+i), where = 1 if site i is occupied with a particle of kind 1, 
Tj = 2 if it is occupied with a particle of kind 2, and Tj = if site i is empty. There are 
M particles of kind 1, and only one particle of kind 2 on the ring. We will refer to them 
as ordinary particles and impurity respectively. 

The system evolves under stochastic dynamics. The possible exchanges between two ad- 
jacent sites during a time interval dt is as follows: 
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The space of configuration is connected. Each configuration can evolve into any 
other and therefore, it has a unique stationary state [18]. As we mentioned above, we 
will perform all the calculations in the grand canonical ensemble and let the number of 
ordinary particles fluctuate around a mean value. We will show that in the large L limit 
the fluctuations in the density of ordinary particles drop to zero so that in this limit, the 
results of the canonical ensemble agree with those of the grand canonical ensemble. 
The above process can be mapped onto the one introduced in [12] by interchanging the 
impurity and vacancies ( 2 <^=^ 0). It can also be considered as a special case of the 
three-species diffusion problems introduced in [13]. 

3 The Stationary Measure and the Quadratic Alge- 
bra 

According to the matrix product formalism, the stationary probability P({r}) of any 
conflguration {r} can be written as a trace of a product of non-commuting operators. 
Since this model is translationally invariant, one can always keep the single impurity at 
site L+1 and write the normalized weight P({t}) in terms of three operators D, E and A: 

P{{r}) - P{{ri,T2,...,TL,TL+i = 2}) = ^Tr[f[{xTiD + y{l-n)E)A] (1) 

The non-zero real variables, x and y, are the fugacities of ordinary particles and va- 
cancies respectively. In the grand canonical ensemble approach, one must choose the 
value of the fugacities to flx the density of ordinary particles to be p = Although it 
would be sufficient to introduce only one fugacity ( since the total number of sites " L+1" 
is fixed), we use both of them to make the symmetry between ordinary particles and va- 
cancies more apparent. The normalization factor Zl in the denominator of Eq.(l), which 
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plays a role analogus to the partition function in equilibrium statistical mechanics, is a 
fundumental quantity and can be calculated using the fact Y^{t} Tl, Tl+i}) = 1- 

Thus one finds 



Zl= E Tr[f[{xnD + y{l-T,)E)A]^Tr{C'^A) (2) 

all configurations i=l 

in which C — xD + yE. The operators E and A satisfy the following quadratic 
algebra: 



DE-qED = D + E (3) 
(3DA = A (4) 
aAE = A (5) 

Assuming that A — \V){W\, Eqs.(3-5) can be written in terms of \V) and {W\ as 
follows 

DE-qED ^ D + E (6) 

^1^) = ^\V) (7) 

(W\E = -(W\ (8) 
a 

With this assumption, the calculation of Tr{...) in expression (2) reduces to the 
calculation of a matrix element 

Zl = (WlC^lV) (9) 
In the next section we will investigate all possible representations of (6-8). 

4 Representations of the Quadratic Algebra 

In [14] the Fock representation of the most general quadratic algebra has been obtained. 
It has been shown that in order to have a finite dimensional representation of certain 
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quadratic algebras, the parameters of the model, say a ,(3 and q in (6-8), should satisfy 
a set of constraints. 

The algebra (6-8) has a one-dimensional representation for a,(3,q < 1, where the oper- 
ators D, E and A are represented by real numbers. Here the vectors \V) and {W\ can 
be discarded and one is dealing with a scalar product state in (9). This representation 
exists if 
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(10) 



and the following constraint holds 



a+P+q^ 1 



(11) 



In Fig.l we have ploted (11) in the three-dimensional parameters space. One can easily 
check that the following matrices [15] 
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and the vectors 
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(14) 



in which 



a — 



(15) 



constitute an n-dimensional representation of (6-8) for a, I3,q < 1, provided that the 
following constraint between the parameters a, {3 and q holds 



a 



l-n 



(16) 



(3 ' ' a 

Note that for n = 1 the above constraint reduces to (11). Using (16) it can be verified 
that for g < 1, the region of the phase space which is accessible by the totally of all finite 
dimensional (n > 2) representations is 

a + (3 + q<l (17) 

In Fig. 2 we have plotted (16) for two values of n. One can see that as the dimension 
of representation increases from 2, the two-dimensional surface (16) begins to approach 
towards the q axis. Therefore the finite dimensional representation of the algebra (6-8) 
allows us to derive exact results for this model only on some special surfaces given by 

(16) . However using the same conjecture proposed in [15], one can introduce a kind of 
analytical continuation to obtain results valid in the whole accessible region given by 

(17) . 

In what follows we show that in the region specified by (17), only three phases can 
exist and this result is independent of n (dimension of the representation). We now find 
the asymptotic behavior of (9) for large values of L. Note that the fugacity of ordinary 
particles x and vacancies y remain constant in this limit. In the thermodynamic limit 
(L — )• oo), the behavior of Zl is governed by the largest eigenvalue of C. The eigenvalues 
6 of 
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C ^xD + yE^ 
f {x + y)+xa + l 



l-q 

can readily be computed 



V 







(18) 



y + y) + xag" + j 



n 



(19) 



ii^-^{{x + y) + xaq'-^ + -^} , z = 1, 

We notice that all the eigenvalues of C lie on the curve [15] 

1 ,, , y. 

z — > [[x + y) + xz + -) 

1 — q z 

in which z — a,aq, aq'^~^. Therefore apart from the values of the fugacities, the largest 
eigenvalue ^^ax of C takes one of the following values: 

/) Cma. = 6, if 6 > Or cquivalcntly x(i^) > yC-^) 

II) Ua. = Cn, if Cn > 6 01 cquivalcntly xC-f^) < vC-^) 

III) Cma. = 6 = Cn, if 6 = Or cquivalcntly x(i^) = y(^) 

These three different cases correspond to different phases. The process is invariant 
when the direction of motion is reversed and the following transformations are applied 



Particle of kind 1 — > Particle of kind (Vacancy) 
Site number i — > Site number L + 1 — i 

P — ' 1 -P 
a — > (3 

Now considering that the interchange of x and y is equivalent to the exchange of the 
density of ordinary particles and vacancies, (/) is symmetric to (II) under these trans- 



7 



formations. 

By using (18) one can see that in the regions (/) and (//), C is diagonahzable, so 
we expect that all correlation functions of form {Ti^...Ti^^-^) depend exponentially on the 
distances involved. In the region (///) the eigenvalues of C coincide 

- Cn-.+i - + y) + y( ^ ~ ^ ~ (/-^ + g"-^)} , = (20) 

1 — g a 

In this case, since C has an off-diagonal line (see 18), it is not diagonahzable. Hence it 
implies an algebraic behavior of all correlation functions (see discussions in [16]). In the 
next section we will show this explicitly for a two-dimensional representation. 
The algebra (6-8) also has infinite dimensional representations [13] and by using them 
one has access to the entire phase space without any constraints on the parameters. Since 
for these representations the calculations seem to be very difficult, we shall adopt one 
and two-dimensional representations to study the general behavior of some interesting 
quantities, such as the density profile of ordinary particles and the speed of impurity in 
different phases. 



5 One and Two Dimensional Representations. Exact 
Results 

As mentioned in the previous section, using n-dimensional representations (n > 2), only 
the region a-\-j3-\-q< \ is accessible, where three phases exist. One-dimensional repre- 
sentation (10) limits us to that part of phase space given by the constraint (11). On this 
two dimensional surface, the partition function (9) has a simple form 

^L = (| + ^)^ (21) 
p a 

The mean value of density of the ordinary particles which is defined by 



8 



L ax 



(22) 



can easily be evaluated 

P - Jrv (23) 

If Pi denotes the expectation that site i is occupied with an ordinary particle in the sta- 
tionary state, knowing that the only impurity is at site L + 1, we have 



Tr{C'-^DC^-'A) _ {W\C'-^DC^-'\V) 
~ ^ Tr{C^A) ~ ^ {W\C^\V) ^ ^ 

Consequently, using (10) and (23), we are able to evaluate the density profile of ordinary 
particles in this case 

Pi^P (25) 

In Fig. 3 we have plotted the phase diagram of the model for a constant q. The accessible 
region (17) lies under the line a + fi = 1 — q. On this line each configuration has the same 
stationary probability and mean field results become exact. In the stationary state, the 
speed of impurity v can be expressed as 

^ yJliEC^ _ Tr{C--'DA) ^ . _ .36) 
^ rr(C^) ^ rr(C^) ' Zl ^ ' 

when use has been made of the algebraic relations (6-8). Using (21), (23) and (26) the 
speed of impurity is found to be 

v = a-{l-q)p (27) 

One may also note that for a > {1 — q)p, (27) becomes negative i.e. the impurity 
starts moving backward. This can easily be understood from the fact that whenever an 
ordinary particle appears behind the impurity, it overtakes the impurity while pushing it 
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backward (Sec section 2). Thus, when the density of ordinary particles exceeds a critical 
value pc = Yz^, a macroscopic negative current of impurity evolves in the system. 



Now we will consider a two-dimensional representation. The following matrices 

1 . \ / 1 ^ 

a 

i a 

'W J \ a 

with the vectors 



D=[ ^ ] , E={ ] (28) 



satisfy (6-8), if the following relation holds 

{q + a){q + /3)^q (30) 

For n — 2, this representation is equivalent to (12-14) and (30) concides with (16). 
As we saw, the properties of the matrix 



i. _|_ 1!. .21 

is of prime importance in determining the phase structre of the model. The eigenvalues 
of C can easily be computed 



6 = | + + , ^ = (1 + 1)^ + ^ 

pa pa 



These eigenvalues also coincide with those obtained from (19) for n — 2. For 7^ ^2, 
one can introduce 

1 / ^ l-g-g 

J- ./n y 



'^(y(i^)-a;(i^l - 
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to diagonalize C 

UCU 



-1 Ki 



V 6 ; 

and using (9) and (29) obtain the following expresion for 



Zl - __fl-q-l. _ „Jl-'^-a^ (33) 



In the thermodynamic limit, two cases can be distinguished. The first case, corre- 
sponding to the phase I, is specified with > ^2, where the first term in (33) becomes 
dominant. The mean density of ordinary particles (22) can be evaluated as 



X 



The inequality ^1 > ^2 can be written in terms of p as 

a> (1-q)(1-p) 

which with a + (3 + q < 1 define the boundaries of the phase (/) (See Fig.3). In the q ^ 
hmit, this phase corresponds to the phase (VI) of the model studied in [1]. Using (24), 
(28-30), (32) and (34) we obtain 

pi = 1 ^(1-p) , Pl = p (35) 

q + a 

The density profile decreases exponentially from pi to p^. One can define a charac- 
teristic length measuring the range of the effect of the impurity 



Note that the correlation length diverges as p — > 1 — ^3^- For L » ^ the density 
profile of the ordinary particles is of the form 
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Pi = cie s + C2e « + C3 (37) 



in which 



Cl = 



C2 



(^1^ - a:i^)6 
a;(?/ — x) 

^ ^(i/^-^i-f^)ei 

Using (26), (33) and (34) we find the following expression for the speed of impurity 
in phase (/) 

" ;i - q)p (38) 



q + a 



One observes that for p > the average speed of the single impurity again 

becomes negative. It can easily be checked that in g ^ limit the values of pi, pi and 
V approach to their corresponding values in the phase (VI) of [1]. Next we will examine 
the second phase which is characterized with ^2 > ^i- This phase is symmetric to phase 
(/) and all the results can be obtained using the transformations introduced in section 
4. Finally the third phase (///) occurs when = ^2, where C is not diagonalizable (See 
32). Nevertheless, Zl can be computed noting that in this case C can be written as 

xi 

xi 




f,-.. 1 - (yi-^i'-^))'^^ (39) 



in which 



5' = (40) 
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Using this property with (9) and (29) yield 

= (W\C-\V) = iA'-^) + ^)-'{-(^) + 1 + LA^^]) (41) 

pa pa p 

and the mean value of density of ordinary particles (22) is found to be 

^_ Vgx+i(l + g)yxy ^^^^ 



+ + (1 + 

The boundaries of the phase (///) (See Fig. 3) are hmited to 

a < (l-p)(l-g) 
/3 < p{l-q) 



which can be distinguished for g < 1 using the fact that x{ ^~^~^ ) = y{ (30) and 

(42). The density profile of ordinary particles, as mentioned above, has an algebraic 
behavior in this phase and increases linearly from pi = to pi — 1 — according to 

p(z) = -J- + {l- , 0<z<l (43) 

For q = this phase corresponds to phase {V) in [1] where a shock structure has been 
observed for certain values of the density p. The linear profile which is the consequence of 
a fluctuating shock front, has also been observed in ASEP and PASEP with open bound- 
aries [11,14,15]. Here the relation x{^—^^) — j/( ^~^~" ) prevents fixing the fugacities x 
and y and we cannot adjust the density p given by (42) to the desired vahic. Therefore 
we cannot see a real shock structure in the grand canonical ensemble. Nevertheless one 
can observe it in this phase using the canonical ensemble [19]. The speed of impurity 
(26) in this phase is of the form 

v^a-(3 (44) 



which is independent of p. Note that for (3 > a the speed of impurity again becomes 
negative. One can check that in this phase the values of pi, pl and v approach their 
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corresponding values in the phase (V) in g — > hmit. 

The fluctuation in the density of ordinary particles [17] which is given by 

can also be calculated in each phase. Using (33) and (41) one can easily check that in 
large L limit these fluctuations drop to zero like ^ therefore the results of the grand 
canonical ensemble agree with those obtained from the canonical one. 



6 Comparison and Concluding Remarks 

As we mentioned in section 2, this model can be considered as a simple generalization 
of the model studied in [1]. Here we compare the results obtained there (we call it the 
model A) with thoes in our model (the model B). 

For q = 0, the quadratic algebra of model B reduces to the one of model A. Using an 
inflnite dimensional representation of this algebra, the author has solved model A exactly 
and shown that it possesses six distinctive phases (I- VI). Since we have used a flnite di- 
mensional representation of the quadratic algebra (6-8), only the region a + f3 + q < 1 lies 
in the accessible region where three phases exist. These phases (I-III) correspond to the 
phases (VI) , (IV) and (V) in model A respectively. In model A, the density proflle of the 
ordinary particles has an exponential behaviour in the phases (VI) and (VI) but it has a 
shock structure in phase (V). In phases (I) and (II) of model B, the density proflle has 
also an exponential behaviour but with a modifled correlation length. In the region (III) 
we have obtained a linear proflle for the density of the ordinary particles in the grand 
canonical ensemble which indicates the presence of a shock in this phase. From the phase 
structure point of view, the backward hopping of the ordinary particles does not change 
the number of phases at least in the accessible region and only the co-existance lines are 
shifted. 

In this paper we have used a flnite dimensional representation of the algebra and done all 
the calculations in grand canonical ensemble which make the calculations rather simple. 
It would be interesting to study the region a + (3 + q > 1, especially the shock structure 
in the phase (III), using an inflnit-dimensional representation of (6-8) in the canonical 
ensemble. 
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Figure Captions 

Fig. 1 Plot of the equation (11). 

Fig. 2 Plot of the equation (16) for n=2 (Left) and n=5 (Right). 
Fig. 3 The phase diagram. 
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Figure 1 . Plot of the equation 
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Figure 3. The phase diagram 



